%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Refraction index of air as a function of wavelength, temperature in degree C, and air pressure in mbar [Edlen].
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;

c0 =299794258;
lmb=[400:1000]*1e-3;														%wavelength

T  =15;																		%temperature in C
P  =1030;                           										%air pressure in mbar

ns=1+0.05792105 ./(238.0185-1./lmb.^2)+0.00167917 ./(57.362-1./lmb.^2);
ns=sqrt(1+4.334446e-4*lmb.^2 ./(lmb.^2-3.470339e-3)+1.118728e-4*lmb.^2./(lmb.^2-1.394001e-2));
nf=1+(ns-1)*0.00185097*P/(1+.003661*T);
 
figure(1);				clf;
subplot(221),plot(lmb*1000,ns-1,'b','LineWidth',1.4)
set(gca,'FontSize',16,'LineWidth',1.4)
xlabel('\fontsize{16}{\it\lambda}   (nm)')
ylabel('\fontsize{16}{\itn}_{r}-1')
grid